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ABSTRACT 

Multi-wavelength observations of the black-widow binary system 2FGL J2339.6—0532 are reported. 

The Fermi gamma-ray source 2FGL J2339.6-0532 was recently categorized as a black widow in which 
a recycled millisecond pulsar (MSP) is evaporating up the companion star with its powerful pulsar 
wind. Our optical observations show clear sinusoidal light curves due to the asymmetric temperature 
distribution of the companion star. Assuming a simple geometry, we constrained the range of the 
inclination angle of the binary system to 52° < i < 59°, which enables us to discuss the interaction 
between the pulsar wind and the companion in detail. The X-ray spectrum consists of two components: 
a soft, steady component that seems to originate from the surface of the MSP, and a hard variable 
component from the wind-termination shock near the companion star. The measured X-ray luminosity 
is comparable to the bolometric luminosity of the companion, meaning that the heating efficiency is 
less than 0.5. In the companion orbit, 10^^ cm from the pulsar, the pulsar wind is already in particle 
dominant-stage, with a magnetization parameter of a < 0.1. In addition, we precisely investigated 
the time variations of the X-ray periodograms and detected a weakening of orbital modulation. The 
observed phenomenon may be related to an unstable pulsar-wind activity or a weak mass accretion, 
both of which can result in the temporal extinction of radio-pulse. 

Keywords: gamma rays: stars — pulsars: general — pulsars: individual (2FGL J2339.6—0532) — 
X-rays: binaries 
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1. INTRODUCTION 

A milli-second pulsar (MSP) is believed to evolve from 
a dead pulsar in a binary system, via mass accretion 
from a companio n star that has sp un up the pulsar for 
billions of years (jAlpar et alJll98^ . Indeed most MSPs 
are discovered in binary systems, however isolated MSPs 
also exist. The missing link between the isolated MSPs 
and binary MSPs is black widow pulsars in which a suffi¬ 
ciently spun up pulsar is evaporating its companion star 
with its powerful pulsar wind. The prototype of the 
black widow is PSR B 1957+2Q with a 1.61 ms radio pulse 
(jFruchter et al.l[l988l ). So far less than 10 black widow 
class objects have been discovered however the evolution 
process from an ac creting MSP to a r otation-powered 
MSP is still unclear (jRobertsI 1^0111 l2Q13l ). 

Moreover a black widow may be a good probe for deep 
investigation of pulsar wind. Currently young rotation 
powered pulsars are believed to have pulsar wind neb¬ 
ulae (PWNe) made up of relativistic electron-positron 
plasma, as seen in the Crab nebula (jWeisskopf et al.l 
I2QQQI: iHester et al.l l2QQ2f) , however the physical mecha¬ 
nism that generates the pulsar wind is not yet under¬ 
stood. Although we do not know whether the descrip¬ 
tions of the Crab pulsar can be applied to a recycled MSP 
with magnetic field four orders of magnitude weaker than 
that of typical radio pulsars, if this is the case, the black 
widows can provide a crucial chance for us to probe pul¬ 
sar wind via the direct interaction with the companion 
nearby the light cylinder. Theoretically pulsar wind is 
dominated by Poynting flux at the light cylinder, there¬ 
fore the magnetic energy in the wind must be converted 
into kinetic energy just after the wind flies out from 
the light cylinder as reported bv lAharonian et ahl (j2Q12f ). 
and we can approach the origin of the pulsar wind much 
deeper with black widows for investigating the unknown 
physical mechanism which can explain the a paradox. In 
order to clarify the history of MSP formation and also 
to constrain the physical mechanism of particle acceler¬ 
ation just around the pulsar, a black widow pulsar is an 
intriguing target. 

The Large Area Telescope (LAT) on the Fermi satel¬ 
lite, with unprecedented sensitivity and angular resolu¬ 
tion in the energy range of 100 MeV to greater than 300 
GeV, has discovered more than 2000 gamma-ray sources 
since its launch, and 30% of these sources are still uniden¬ 
tified. A bright gamma-ray source discovered at high 
galactic latitude vicinity, 2FGL J2339.6—0532 (IFGL J 
2339.7—0531), was also listed in the first Fermi source 
catalog as an un-iden tified source (|Abdo et al.l 120101: 
lAckermann et al1l2012f) . The gamma-ray flux amounts 
to 3.0 ± 0.2 X 10“^^ erg s“^ cm“^ with a variable in¬ 
dex of 15.7, indicating that the gamma-ray flux seems to 
be steady at the month time scale. While the gamma- 
ray spectrum has a cutoff structure at 3 GeV. These 
characteristics lead us to believe the object is a pulsar. 
A follow-up X-ray observation conducted with Chandra 
discovered an X-ray point source wit hin the error circl e 
expected from the gamma-ray image (jKong et ani2012f ). 
However a radio pulse was not detected at the position 
at that moment. On the other hand, ground based op¬ 
tical observations discovered clear sinusoidal variability 
with a period of 4.63 hours, which implies the object 
is in a binary system and the observed optical variabil¬ 


ity is likely related to the orbital motion. The intensity 
of the optical counterpart varies from 20 to 17 mag in 
R-band (jRomani fc Shawl[2Qlll:iKong et al.ll2Q12[ ). More¬ 
over the phase-resolved spectroscopy indicated that the 
companion may be K-class star with a mass of 0.075 Mq. 
This means that the pulsar side hemisphere of the com¬ 
panion star is drastically heated and the energy might 
be supplied from an unknown recycled MSP via pulsar 
wind, as seen in PSR B1957-1-20. 

In this paper, we report multi-wavelength observations 
of the newly discovered black widow binary system 2FGL 
J2339.6-0532, covering near-infrared, optical and X-ray 
energy band. In section 2 and 3, the optical and X- 
ray observations and the obtained results are described, 
respectively. In section 4, we discuss the orbital param¬ 
eters based on the obtained phase resolved SED and the 
properties of the pulsar wind just around the pulsar as 
well as an interpretation of the X-ray light curve showing 
intriguing irregular variability. 

2. OPTICAL OBSERVATIONS 

2.1. Observation 

2FGL J2339.6—0532 was observed from September 22 
to October 7 2011 utilizing the global telescope network, 
“Optical and Infrared Synergetic Telescopes for Educa¬ 
tion and Research (OISTER)” E3- OISTER consists of 14 
independent observatories that are funded by Japanese 
universities and research associations. Eor this work, we 
also asked for photometry observation from the Kottamia 
Astronomical observatory which has a 188 cm reflector. 
Thanks to the locations of the observatories distributed 
across the globe, we can continuously trace the variabil¬ 
ity all day long. Moreover the telescopes cover a wide 
range of wavelengths from Ks to B band which can pro¬ 
vides phase resolved spectral energy distribution (SED) 
of the target. The conducted photometric observations 
are summarized in Table IH The locations of the reference 
stars and the target are shown in Eigure[TJ 

2.2. Ha imaging 

In the case of PSR B1957+20, the binary system is 
surrounded by an Ha bow-shock due to its super sonic 
proper motion with res pect to the surrounding ISM 
(jKulkarni fc Hestei1ll988f ) . To search for a similar feature 
around 2EGL J2339.6-0532, we took Ha images with the 
Kiso Schmidt telescope at the minimum phase to prevent 
contamination from the companion star. Eigure[T] shows 
the obtained Ha image in the vicinity of 2EGL J2339.6- 
0532 utilizing an Ha narrow band filter (AA = 9 nm), 
in which we could not find the evidence of bow-shock 
nebula. 

To evaluate the detection limit, we calculated the stan¬ 
dard deviation of the 33"x67" (50x100 pixel) rectangu¬ 
lar region near the target and obtained ST DEV = 7.82 
ADU. The sky background of the GCD image shows 
a weak striped pattern running in the southeast-to- 
northwest direction, with a peak-to-peak value of ^ 8 
ADU, which seems to be the dominant noise compo¬ 
nent restricting the detection limit. Regardless, if we 
adopt a conservative threshold with a 3cr confidence 
level, the detection limit becomes 23.5 ADU. Based on 

http://oister.oao.nao.ac.jp/ 
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Figure 1. Ho; image of the vicinity of 2FGL J2339.6—0532 observed with the Kiso-105cm Schmidt telescope. The circle region with a 
radius of lO^indicates the location of the target. Two rectangular regions indicate the photometric reference stars shown in Table [T] 


the flux calibrations of the reference stars listed in Ta¬ 
ble [U the conversion coefficient of ADU to energy flux, 
(8.3 ± 0.1) X 10“^^ erg s“^ cm“^ ADU“^, was observed 
for the Ha filter with a bandwidth of 9 nm. Finally, we 
obtained a 3cr detection limit of < 1.9 X 10 erg s ^ 
cm“^ pixel“^, which corresponds to a surface brightness 
of < 8.7 X 10“^^ erg s“^ cm“^ arcsec“^ E3- 

The prototype black widow PSR B1957+20, residing 
at a distance of 1.6 kpc from the earth, has a bow-shock 
nebula extending over ^ 30" x 30"with an Ha flux of 
3.3 X 10erg s~| inte grated over the entire nebula 
(jKnlkarni fc Hest^ll988F If our target is accompanied 
by an equivalent bow-shock nebula, the Ha flux and the 
size should be 6.1 x 10“^^ erg s“^ cm“^ and ^ 41" x 41", 
respectively. Therefore, the expected surface brightness 
is 3.6 X 10“^^ erg s“^ cm“^ arcsec“^, which is slightly 
lower than the detection limit. 

Compared with the other bow-shock PWNe, the size 
and intensity both seem to be scattered across a range of 
two orders of magnitude, possibly refle cting their shock 
conditions (jChatterjee fc CordesI l200^ . Therefore, we 
conclude that this observation provides only a weak up¬ 
per limit of surface brightness in the Ha band. 

2.3. Photometry 

The other telescopes except for the Kiso observa¬ 
tory carried out photometric observations covering Ks^B 
bands from September 9th 2011 to September 30th 
2011. For flux calibration we chose two reference stars 
near the target object, U0825_19993817 (R=17.5) and 
U0825_19993871 (R=16.4) in USNO-2.0A catalog, so 
that all of the observatories can employ the common ref¬ 
erences without saturation. The photon fluxes of the 

Kiso 2k CCD’s pixel resolution is 1.5^^ (IItoh et al.|[200Th . 


references stars were measured and co mpared to the 
standard stars in the Landolt catalog (|Landoltl Il992[) 
with MSI of the Pirka telescope in Hokkaido prefecture 
(Japan). Since the Pirka employs the Johnson-Cousins 
filter system, we estimated the photon flux in the SDSS 
g’-band based on the observed B an d V magnitudes usi ng 
a conversion equation proposed by iSmith et al.l (j2QQ2[ ) , 

g' = VP 0.54(H -V)- 0.07. (1) 

For J, H, and Ks, we employed the photometric data in 
the Two Mi llion All Sky Survey (2 MASS) catalog for flux 
calibration (jSkrutskie et al.l[2006F The obtained magni¬ 
tude of the reference stars are summarized in Table [TJ 

Figure [2] shows the light curves of 2FGL J2339.6—0532 
obtained by OISTER. We re-co nfirmed the clear sinu¬ 
soidal modulatio n as re ported by iKong et ahl (j2012[ ) and 
iRumani fc Shawl (j201l[ ) . The absolute flux seems consis¬ 
tent with the past observations: the modulation ampli¬ 
tude is about 4.5 magnitude in the R-band and is larger 
at shorter wavelength. In R-band, we successfully ob¬ 
served the maximum phases of the target 8 times dur¬ 
ing the observation campaign. The photon flux at the 
maximum phase is Rmax = 17.7 mag on average and 
varies among peaks within a range of ±0.1 mag. There 
was no irregula r activity like the Har es observed in 2FGL 
J1311.6-3429 (|Kat,a.oka, et, a1.ll2nT^ . 

2.4. Phase-resolved SED 

To discuss the energetics quantitatively, we estimated 
the energy flux based on the photometric data shown 
in Figure [21 First we correct ed the galactic extinction at 
the coordinate of th e target (jSchlafiv fc Finkbeinerll20lTI: 
iCardelli et al.llTosM^ : the employed extinction at vari- 

In the above extinction evaluation we utilized 
http://ned.ipac.caltech.edu/forms/calculator.html, and 
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Table 1 

Summary of Field photometry 



Reference-1 

Reference-2 

AAt 

Name(USNO-2.0A) 

Coordinate(J2000.0) 

U0825-19993817 
(23:39:39.487, -05:32:40.56) 

U0825_19993871 
(23:39:40.366, -05:31:51.89) 


B 

18.435 ±0.014 

17.603 ± 0.008 

0.120 

V 

17.834 ±0.018 

16.828 ± 0.008 

0.091 

R 

17.484 ± 0.016 

16.388 ± 0.007 

0.072 

I 

17.296 ± 0.020 

16.105 ±0.008 

0.050 

g’ 

18.098 ±0.011^ 

17.177 ±0.006^ 

0.110 

J 

16.47 ±0.10^ 

15.19 ±0.04^ 

0.024 

H 

16.45 ±0.23^ 

14.87 ±0.05^ 

0.015 

Ks 

16.12^’^ 

14.55 ± 0.08^ 

0.011 

Ho; (656 nm) 

314 ± 46 

854 ± 170 



Note. — Errors are with Icr confidence level. 

^ Galactic extinction at the target coordinate dSchlafiv fc Finkbeinerll2Qlll: ICardelli et al.iri989f) 

^ SPS S magnitudes were calculated from B-band and V-band magnitudes based on ISmith et al.l 
(l2QQ^ . 

^ For the fiux calibrations of J, H, and Ks, we referred the 2MASS catalog dSkrutskie et al]|2Q06r) . 
^ The catalog error of the Reference-1 for Ks-band was not available. 

^ The fiux density was evaluated from the absorption corrected SED assuming a simple balckbody 
model. The estimated temperature of the reference stars were TBB,refi = 7320 ± 160K, TBB,ref 2 = 
6300 ± 180 K. 


ous wavelength are listed in Table [H Then we converted 
the absorption corrected magnitu des to energy flux using 
con version equations presente d in iFukugita et al.l (jl996h 
and iTokunaga fc Vaccal ()2QQ5f ) for optical (B, V, R, I, g’) 
and IR(J, H, Ks) energy bands, respectively. 

Figure [3] — (a) shows the yielded energy fluxes as func¬ 
tions of the binary orbital phase. Phase=0 was set to 
be MJD=55500.0. While the optical light-curves show 
symmetric structures, IR (H and Ks) light-curves pos¬ 
sess somewhat asymmetric shapes. This may be due to 
the geometry of the companion star or evaporating stel- 
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Figure 2. Multi-color light curve of 2FGL J2339.6—0532 obtained 
via OISTER. 


http://dogwood.physics.mcmaster.ca/Acurve.html 


lar gas. Panel (b) shows the spectral energy distribution 
during the brightening phase from orbital phase= 0 to 
0.5 in panel (a). Clearly the peak frequency increases as 
the orbital phase increases from orbital phase 0 to 0.5, 
implying that the temperature of the companion is in¬ 
creasing with the orbital phase. At the maximum phase 
the SED is fitted with a black-body model with an ef¬ 
fective temperature of Teff = 7540 T 130 K, which is 
consistent with the past study (jRGmani fc Shawll201l[ ). 

In this paper, we assumed a simple emission model 
from the surface of the companion star to explain the 
observed optical emission. As described in Figure IH we 
supposed that a companion star with a radius Tcomp is 
orbiting around a pulsar at an orbital radius Rorb, and 
the hemisphere of the companion that faces the pulsar 
is heated by the pulsar wind, while the opposite side of 
the companion has ordinal temperature Tcooi- Note that 
this model does not take into account the effect of energy 
transfer via convection or advection on the surface of the 
companion, therefore the temperature distribution of the 
companion is simply described by the energy injection via 
pulsar wind although the heating efficiency is unclear. In 
this paper, we assumed that the spin-down energy is per¬ 
fectly converted into the isotropic pulsar wind and the 
injected energy via the pulsar wind heats the compan¬ 
ion up with a heating efficiency of / for simplicity. We 
also adopted a spin-down luminosity of 2.3 x 10^^ erg 
s“^ based o n the timing ana lysis in radio band recently 
reported by iRnv et all (j2Q14f ) . 

We then fitted the obtained phase-resolved SED with 
the model function described above (Eigure[5]). Since the 
model seems to correlate weakly with the inclination an¬ 
gle, we fixed the inclination angle i. The resultant model 
parameters are summarized in Eigure [6] as functions of 
inclination angle. Although the fit seems rather poor — 
mainly because of the large residual in the IR band — we 
found that an inclination angle i = 59° minimizes the 
Eor comparison, we also attempted this analysis without 
IR data and obtained a best-fit parameter of i = 52°. 
These resul ts seem consisten^ t with the inclination angle 
reported bv iRGmani fc Sha^ (|2QllF The obtained best- 
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(a) Phase folded light-curves 


(b) Spectral energy distributions at increasing phase 



I phase = 0.4^.5 

: phase = 0.3~0.4 









piiciac = u.z 

- phase = 0.1 

: phase = 0.C 

l~0.2 

)~0.1 









- 




— 





- 









: - 

-i- 

. .4^. 







U 

I I I I 


T 



,~f 

i * *i 


i i 



0.8 1 1.2 1.4 

Orbital Phase 


3 4 5 

Frequency [10^^ Hz] 


7 8 9 10 


Figure 3. (a)Phase folded multi-color light curves. The assumed orbital period is 4.63435 hr and Phase=0 is set to be MJD=55500 day. 
Marker colors represent the energy bands: red= R, green= V, blue=B, magenta=H, and black=Ks. (b)Spectral energy distributions at 
increasing phase from the minimum to the maximum. Colors correspond to the orbital phase: black=0-0.1, magenta=0.1-0.2, red=0.2-0.3, 
green=0.3-0.4, and blue=0.4-0.5. 



Figure 4. Schematic image of the assumed emission model in 
which a companion star with a radius rcomp is orbiting around a 
pulsar at an orbital radius Rorb • The hemisphere of the companion 
star facing the pulsar is heated. In this paper, the inclination angle 
i is dehned as the angle between the line of sight and the normal 
line to the orbital plane. 

fit parameters are summarized in Tabled Additionally, 
the obtained temperature distributions of the compan¬ 
ion surface are plotted in Figure 0 as functions of the 
zenith angle of the pulsar from the measured point on 
the companion surface. 

3. X-RAY OBSERVATION 
3.1. Observation 

2FGL J2339.6—05312 was observed with the Suzaku 
satellite on Jun. 29th, 2011. The target was observed at 
the XIS nominal point with the normal imaging mode. 
The net exposure time after the standard data reduction 
process was 96 ks. The accumulated photons are 2555 
counts with the front-illuminated CCDs (XISO + XIS3) 
and 1874 counts with the back-illuminated CCD (XISl) 
in the 0.4—8 keV energy band. During the observation, 
the event rates were ^ 8.25 ± 0.28 x 10“^ (XISO, XIS3) 
photons s“^ and 8.89 ± 0.50 x 10“^ photons s“^ (XISl) 
in the energy band. 

3.2. Spectral analysis 

Figure [8] shows the X-ray spectrum obtained with the 
Suzaku. At first, we tried to fit the X-ray spectra with a 



Figure 5. Phase-resolved SED fitted with a model function. 

single power-law function however the resultant fit was 
poor (x^/u = 67.71/55) due to residuals below 1 keV, im¬ 
plying the presence of a soft excess component. We there¬ 
fore added a black-body to the model function for the 
soft excess, and obtained a better fit (x^/^ = 51.82/55). 
Dotted lines on Figure [8] describe the best-fit model func¬ 
tions. The resultant parameters of the model fitting 
are summarized in Table [3l In both case, the column 
density converged on zero. The resultant upper limit is 
< 1.2 X 10^^ cm“^ with 90 % confidence level, and is con¬ 
sistent with th e Galactic absorption of 3.2 3 x 10^^ cm^ for 
this direction (jDickev fc LockmanI 119901 iKalberla et al.l 

inni). 

In order to clarify the origin of these two components, 
phase-folded X-ray light curves were generated as shown 
in Figure [3 In this figure, Phase=0 was set to be 
MJD=55500.0 in the same way as the optical light curves 
shown in Figure [3l-( a). Panel (a)^(c) in Figure [9] rep¬ 
resent the energy bands 0.4—1.0 keV, 2.0—4.0 keV and 
4.0—8.0 keV, respectively. Clearly, the hard X-ray light 
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Table 2 

Obtained parameters of the binary system from the phase-resolved SED 


Data set 

Optical+IR^ 

Optical^ 

Orbital radius (Rorb)^ 

1.08 ±0.02 X lO^i cm 

1.18 ±0.02 X lO^i cm 

Companion radius (rcomp) 

1.26 ±0.02 X lO^o di 1 cm® 

1.53 ±0.02 X lO^o dl 1 cm® 

Companion Temperature (Tcool) 

3170 ±70 K 

2740 ± 120 K 

Heating Efficiency {f)^ 

1.42 ±0.05 

1.16 ±0.04 

Inclination angle (i) 

59.°0±1.°5 

52.°1 ± l.°0 

Bolometric Luminosity® 

1.09 ± 0.17 X 10^2 dl 1 erg s"! 

1.12 ±0.16 X 1032 dl 1 erg s"! 


1392.6 

637.6 


Note. — Errors are with 1 a confidence level. 

^ “Optical” and “IR” correspond to [B, V, R, I, g’] data-set and [J, H, Ks] data-set, respectively. 

^ Rorb are calculated based on radial velocity of the companion star reported bv iRomani Sz Shawl (I2011f) . 

^ di.i is the distance in unit of 1.1 kpc. _ 

Assuming an isotropic pulsar wind with a luminosity of 2.3 X 10^^ erg s~^ based on I Ray et al.l (1201411 . 

® Bolometric luminosity of the companion’s heated hemisphere. 


Table 3 

Results of model fitting of the averaged X-ray spectra 


Model 

wabs*(bbody ± power-law) 

wabs*power-law 

A^h (10^^ cm-2) 

< 1.2 (0) 

< 0.13 (0) 

Blackbody Temperature (keV) 

0.15 ±0.06 


Blackbody Radius (km)"^ 

0.28lg;“ 


Power-law Photon Index 

6 d 
+ 1 

r—1 

1.32 ±0.08 

Power-law Fluxo. 5 _iokeV (10“^3 gj-g cm“^ 

2.50 ±0.16 

2.50 ±0.15 

Total fluxo.5_iokeV (10“^3ej.g cm“2 g-i) 

2.67i;;°i 

73 

(dof) 

51.82 (51) 

67.71(53) 


Note. — Errors are with 90 % confidence level. 

^ Blackbody radius assuming a distance of 1.1 kpc. 

curves show clear modulation coinciding with the orbital 
motion. This timing coincidence may imply that the non- 
thermal emission originates from the companion surface. 
However the shape of the light curves is somewhat differ¬ 
ent from the sinusoidal shapes observed in optical, i.e., 
a double-peaked shape at 2.0—4.0 keV band, and a flat 
top shape at 4.0—8.0 keV. On the other hand, the soft 
X-rays seem to be steady. This may imply that the soft 
component has a different emitting region, namely the 
pulsar. 

To evaluate the upper limit of the modulation fac¬ 
tor in the soft X-ray energy band, we tried to fit the 
X-ray light curve with a sinusoidal curve and a con¬ 
stant component, Hsin(27rx) + H, where x is the orbital 
phase. The obtained amplitude of the sinusoidal curve 
is H = (0.80 ± 0.34) x 10“^ counts s“^ against the con¬ 
stant component of B = (3.95 ±0.23) x 10“^ counts s“^. 
Therefore, the upper limit of the modulation factor is 
AjB < 28.6% with 1 a confidence level. Note that the 
maximum point was not fixed to be phase=0.5 in this es¬ 
timation; this helped to maximize the modulation factor. 
The best-fit model curve is shown in Figure [9] (a). 

If we adopt the two-component model consisting of a 
blackbody and a power-law function from the X-ray spec¬ 
troscopy, we can roughly evaluate the contributions from 
these two components to the constant emission. The 
phase-averaged absorbed-photon flux for the blackbody 
and the power law in the 0.4—1.0 keV energy band are 
both 2.0 X 10“^ photons s“^. Assuming that the con¬ 


stant component originates in the blackbody, the modu¬ 
lation factor is expected to be {PL) j{BB ^ PL) ^ 50 %, 
which is about twice of the measured value. Taking into 
account the large errors arising from the uncertainties 
of the blackbody components (see Table m, the mea¬ 
sured modulation factor should be larger because of the 
tight upper limits to the blackbody component. This 
may indicate that the constant emission may consist of 
the blackbody from the NS surface and the non-thermal 
component originating in the pulsar magnetosphere. 

3.3. Temporal analysis 

The upper panel of Figure [10] shows the obtained X- 
ray light curve with a bin size of 600 s for the energy 
range of 0.4—8.0 keV. Note that the whole length of 
the light curve is 190 ks, containing occultation peri¬ 
ods of earth, and corresponding to the twice of the net 
exposure time of 96 ks. Although the photon statistics 
seem poor, the X-rays show periodic variability coinci¬ 
dent with an orbital period of 16683 s. We expected to 
observe the flare-like f eatures reported in previous studies 


(iKataoka et al.ll2012l: iPa 

Ditto et al.ll2013l:lFerrimo et al.l 

12014 iLinares et al.ll2014 

), which may be related to the 


accretion activity; however, no significant evidence for 
this was discovered. Nonetheless, the shapes of the light 
curve seem to change over time. 


In order to clarify the temporal variation in the wave¬ 
form of the X-ray light curve, we traced the temporal 
changes of the power spectrum with the Lomb-Scargle 
algorithm, which can be applied to an unevenly sampled 
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Figure 6. The obtained parameters from model-fitting of the 
phase-resolved SED as functions of inclination angleL Panel 
(a)~(d) show the orbital radius Rorb and companion radius rcomp, 
the companion temperature at the cool side, the heating efficiency 
/, and the resultant (d.o.f = 901), respectively. The gray-bold 
line on panel (a) represents the semi-major axis of the compan¬ 
ion calculated from the radial velocity of 350 km s~^ with re spect 
to the mass centroid of the system (|Romani &: Sha^ 120111) . In 
this calculation we employe d a spin-down lum inosity of the pulsar, 
Lsd = 2.3 X 10^^, based on IRav et al.l (120141 b 

data set. In this analysis, we employed a rectangular 
function as a window function to detect small changes 
from a short data sequence. The width of window func¬ 
tion was set to 50 ks in order to cover three orbital laps 
and to detect the light-curve variations for the total ex¬ 
posure time of 190 ks. We then scanned the light curve 
by sweeping the window function, and generated a power 
spectrum for every timespan in the light curve. Note that 
this method cannot detect rapid variations of the power 
spectra (faster than a few tens of kilo-seconds), because 
the width of the time window has a finite timespan. 

The lower panel of Figure [10] shows the obtained pe- 
riodogram as a function of time, which represents the 
central time of the sampled time window. Therefore, 
the periodogram at t = 25 ks consists of the X-ray light 
curve in the time span of 0 < t < 50 ks. Both ends of 
the grayscale map are omitted because the time window 
can sweep in the range of the obtained light curve. The 
sweeping step size was set to 600 s. Therefore, the time 
resolution of the running periodogram is about 230 which 



Figure 7. Simulated temperature of the companion surface as a 
function of zenith angle of the pulsar from the position for which 
the calculation is performed. The simulations were performed for 
the best-fit parameters listed in Table [2] 



0.5 1 2 5 10 

Energy [keV] 


Figure 8. X-ray spectra observed by the XIS aboard the Suzaku 
satellite. Black and red data points corresponds to the spectra 
observed with the front illuminated CCDs (XISO and XIS3), and 
with the back illuminated CCD (XISl), respectively. Dotted-lines 
show the best-fit model functions consisting of a low-temperature 
black-body and a power-law function. The model parameters are 
found in Table [3] 

is much greater than the timescale that can be detected 
by this method. Nevertheless, a clear peak structure was 
detected at a frequency of around fi = 0.6 x 10“^ Hz, 
which corresponds to the orbital frequency. On the other 
hand, /2 = 1.2 x 10“^ Hz appears to denote the second 
harmonic of the orbital frequency, which may be related 
to the dip structure of the double-peaked light curve. In¬ 
terestingly, these two components in the power spectra 
decrease sequentially after t = 1.2x10^ s. 

4. DISCUSSION 
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Orbital Phase 

Start Time 15741 14:28:32:096 Stop Time 15743 19:28:01:748 

Figure 9. Phase-folded X-ray light curves for (a)0.4—1.0 keV (b)2.0—4.0 keV and (c)4.0—8.0 keV energy bands. Phase=0 is set to be 
MJD=55500.0 in the same way as the optical light curves. Dotted line in panel (a) is the best-fit sinusoidal curve for evaluating the 
modulation factor. Dashed line in panel (c) shows a model function for synchrotron emission from a thin layer covering the companion 
hemisphere. 


4.1. Constraints on the inclination angle of the binary 

orbit 

The model fitting of the phase-resolved SED resulted 
in inclination angles of i = 59° from the entire data 
set (covering optical -h IR bands), and i = 52° from 
the optical-only data. This discrepancy seems to arise 
from the systematic errors of the reference stars in the 
IR band (J, H, Ks), resulting in an underestimate of the 
target luminosity in the IR band. Therefore, the SED 
of the optical-|-IR data set tended to return higher com¬ 
panion temperatures than the optical-only data. On the 
other hand, it is also difficult to constrain the tempera¬ 
ture with the optical-only data (B, V, g’, R, and I), which 
may then result in an underestimate of the companion 
temperature. Therefore, the true inclination angle must 
be in the range of 52° < i < 59°. Hereafter, we take 
i - 55° ±5°. 

We should also summarize the pul sar parame t ers in 
preparation for the discussion below. iRav et aP (|2Q14f ) 
recently reported the detection of radio and gamma-ray 
pulse emission from the central MSP with a period of 
2.88 ms with a spin-down luminosity of Lsd = 2.3 x 10^^ 
erg s“^. Based on these parameters, the pulsar’s sur¬ 
face magnetic field can be estimated, assuming that the 
pulsar’s rotatio n energy is dissipated by dipole radiation 
(jLongairlll994f ). 

= 3.0 X = 1.9 X 10*14^/^ G, (2) 

where is the pulsar’s moment of inertia in units of 10^^ 
g cm^, which is a typical magnetic field for a standard 
MSP. 


4.2. Condition for mass Accretion 

In the case of neutron star (NS) binary systems, the 
mass accretion pro cess is thought to be controlled by 
the propeller effect (jStella et al.lll986l ). e.i., the pressure 
balance between the ram pressure of accreting material 
Pram and the magnetic pressure of the pulsar pb- Here 
we define rm where pe and Pram are in balance. 



(li^) ( iQi^is-O cm, (3) 

where G is the gravitational constant, r^ is the NS radius, 
is the NS mass, m is the mass accretion rate. Nearby 
the pulsar, charged particles are trapped by the pulsar 
magnetosphere and are co-rot at ing with the pulsar. If 
the accreting gas approaches the neutron star within a 
radius. 



= 1-7x10'*(t:^) ( 2 :^)'^' cm (4) 

at which the centrifugal force and the gravitational at¬ 
traction are equal, the gas will accrete on the neutron 
star (where is the angular velocity of the NS). 




























Multi-wavelength observations oe the black widow pulsar 2FGL J2339.6-0532 


9 




o 

O 


N 


>1 

U 

c 

Q> 

3 

CT 

P 



Figure 10. Upper panel shows the X-ray light curve for energy band of 0.4—8.0 keV with a bin size of 600 s. Lower panel shows the 
running periodogram generated from the X-ray light curve shown in the upper panel with the Lomb-Scargle algorithm. We employed a 
50ks rectangular window function that scans over the entire X-ray light curve of 190 ks. The horizontal line in the lower panel corresponds 
to the center time of the sampled time span of 50 ks. The gray scale shows arbitrary amplitude of modulation. Note that the amplitude 
of 6.3 correspond to 90 % confidence level. 


Utilizing from the pulse profile, the lower limit of 
the accretion rate m that fulfills the condition of Vc > 
can be calculatecQ, 

^ 23/2G'5/3mU^ 

Recently, iRnv et aP (j2Q14f ) reported temporal extinc¬ 
tions of the pulse emission in radio implying that the 
accretion on the MSP is ongoing. If that is the case, 
the mass accretion rate must be higher than Eq. [5] 
and the expected luminosity during the mass accretion 
will amount to 2.4 x 10^^ erg s“^ which is about 5 or¬ 
ders of magnitude higher than the X-ray luminosity at 
the 0.5—10.0 keV energy band. iPapitto et al.l (j2013[ ) re¬ 
ported intriguing activities of IGR J18245—2452 in the 
interme diate stage be t ween the r otation and ac c retion 
power. iLinares et al ] (j2014f ) and iFerrigno et al.l (|2014l f 
reported weak accretion flows (below m < 10^ g s“^) 
in IGR J18245—2452; these exhibit a variety of activi¬ 
ties possibly reflecting accretion conditions. In case of 
Tm ~ Tc, namely the weak-propeller regime, accretion 
flow is only partly inhibited and the X-ray luminosity 
is in the range of 10^^ ^ 10^^ erg s“^. For cases of 
namely the strong propeller regime, the X-ray 
luminosity is in the range of 10^^ ^ 10^^ with a very fast 
and striking variability. Regardless, the X-ray luminos¬ 
ity of 2FGL J2339.6—0532 is still lower than that of IGR 
J18245—2452 in the strong-propeller regime by an order 
of magnitude. 

In addition IGR J18245—2452 exhibited another qui¬ 
escent state in 2002 with an X-ray luminosity of ^ 10^^ 


erg s“^, which is comparable to t he observed X- r ay lu- 
minosity of 2FGL J2339.6—0532. IFerrigno et al.l (j2014f ) 
claimed that reached the light-cylinder radius in this 
regime. This condition would allow for the existence 
of pulsar wind and would also prevent mass accretion. 
Thus, we conclude that the binary system can be in the 
same condition seen in IGR J18245—2452 in 2002 or in 
the late stage of the evolution process without mass ac¬ 
cretion. In the latter case, where the mass accretion 
already ended, the radio-pulse extinction can also be ex¬ 
plained by the unstable activity of the pul sar’s magne¬ 
tosph ere itself, as with “nulling” pulsars (jGajiar et al.l 
I 2 OI 2 L see also references therein.). 

4.3. Origin of the X-ray emission 

As shown in the spectral analysis, the soft X-rays be¬ 
low 1 keV seem to originate from a black-body. The 
obtained radius of 0.28 km for d = 1.1 kpc is unusually 
small for an accretion disc and thus the radiation can be 
interpreted as the X-ray radiation from the neutron star 
surface. A radius of 0.28 km is still small for a neutron 
star, therefore the emitting region could be a hot spot 
around the polar cap area. In order to constrain the ge¬ 
ometry of the polar cap, a high time resolution X-ray 
observation is required. In addition, the small variability 
and the lack of eclipse in X-ray light curve are consistent 
with the estimated inclination angle of i ^ 55°. 

In contrast, the hard X-ray light curves vary with the 
orbital motion. The X-ray luminosity of the non-thermal 
component from the averaged X-ray spectra amounts to 
Apow = 3.6 X 10^^ erg s“^ for a distance of 1.1 kpc at 
0.5—10.0 keV. Therefore the ratio of Tpow against the 
spin-down luminosity is about 0.15 %, which is the typi¬ 
cal value for a standard radio puls ar (jBecker fc Trnemperl 
I 1997 I: iKargaltsev fc Pavlovll200^ . 


In which we assumed a relation 745 oc M^r^. 
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To explain the flat top modulation of the 4—8 keV 
light curve, two possible emitting regions can be pro¬ 
posed. One is an extended emitting region just around 
the pulsar like a halo. In this hypothesis the minimum 
phase at the inferior conjunction can be explained by 
an eclipse of the halo. The observed smooth modula¬ 
tion requires the emitting region to be larger than the 
companion, lO^^cm^, which is still small for a standard 
pulsar wind nebula. In many case the non-thermal X- 
rays from a pulsar wind nebula is from the downstream 
of the termination shock and the radius of the termina¬ 
tion shock is constrained by the pressure balance. In this 
case, the upper limit of the emitter size is the orbital ra¬ 
dius of Rorh ^ 10^^ cm. At this distance however the 
wind pressure is still too strong to be terminated. 

The other hypothesis for the origin of the hard X-rays 
is that of a shocked pulsar wind interacting with the com¬ 
panion star. The companion star orbiting nearby the 
pulsar can partially terminate the strong pulsar wind. 
In this case the flat-top light curve can be explained by 
the Lorentz boost of shocked plasma. As shown in Fig¬ 
ure m we assumed a thin emitting layer covering the 
hemisphere of the companion star. The surface of the 
emitting region is the wind termination shock and the 
injected electron plasma decelerated at the shock to a 
bulk velocity /32 reflecting the magnetization parameter 
upstream of the shock. At the superior conjunction phase 
the shocked plasma is flowing in a direction away from 
us, t herefore the synchr otron luminosity may be weak¬ 
ened (jPelling et al.lll98^ . The dashed line in the bottom 
panel of Figure [9] shows the calculated photon flux for the 
thin layer of the shocked pulsar wind on the surface of 
the companion with a photon index of p = 1.14 and an 
inclination angle of i = 55°. The best-fit bulk velocity of 
the shocked wind is /32 = 0.43 ± 0.15 in units of the light 
speed. If that is the case, the X-ray flux is attenuated 
by a factor of 




i+P 


1-/32 cos(7r/2 + i) 


0.42. 


(6) 


Here we employed a magnetization parameter a 
for the pulsar wind upstream of the shock , a = 
H^/(47rnii4i7i, rUeC^) (|Kennel fc Coronitil I1984J) . where 
Bi is the magnetic field just upstream of the shock, ni 
is the number density of the pulsar wind, Ui is the flow 
velocity defined hy 1 -\- (Fi is the bulk Lorentz 

factor of the wind upstream of the shock), and c is the 
light velocity. Assuming that half of the wind energy is 
carried as Poynting flux (cr = 1), Hi amounts to 


Hi 




1/2 


= 5.5 G. 


(7) 


Furthermore, the pulsar wind is compressed at the shock, 
therefore the downstream magnetic field H 2 is higher 
than Hi. Adopting H 2 > 5.5 G as a lower limit, we can 
constrain the gyration radius of the synchrotron electrons 
emitting 1 keV X-rays, 


I gyro 



( B2 \ 

VlkeV/ ' 

V5.5 G/ 


-1 


(8) 


Termination shock 


Contact Discontinuity 



Shock-heated pulsar wind 


Figure 11. Schematic view of the Hard X-ray emission from the 
surface of the companion. 


where ex is the synchrotron photon energy. The ob¬ 
tained gyration radius is far smaller than the compan¬ 
ion, Tcomp ^ 10^^ cm, therefore the shock heated elec¬ 
trons cannot escape from the downstream immediately. 
In terms of the enegetics, the phase averaged X-ray lumi¬ 
nosity of the power-law component was Hpow,o. 5 -iokeV = 
3.6 X 10^^ erg s“^ from the spectral fitting. Taking into 
account the Lorentz boost effect of 1/0.42, and the ratio 
of the maximum flux against the phase averaged flux of 
^3.0, the intrinsic synchrotron luminosity from the com¬ 
panion becomes L'^ow,o. 5 -iokeV = 1-1 x 10^^ which is com¬ 
parable to the bolometric luminosity of the heated hemi¬ 
sphere of the companion, Hbb ^ 1.1 x 10^^ ergs s“^. This 
brief estimation also indicates that the heating efficiency 
should be smaller than 0.5. This result obviously differ¬ 
ent from the heating efficiency obtained from the phase- 
resolved SED (see Figure [5]) in which an isotropic pulsar 
wind is assumed. This discrepancy may indicate that 
the pulsar wind from the cen tral MSP is not isotropic as 
expec ted in the crab nebula (| Bogovalov fc Khangonlvanl 
[2001 ). The wind pressure at the companion is then 
Pwind = Lsb/^T rRlrh^ ^ 4.7 dyu which may be greater 
than that of the coronal gas (T = 10^ K, Up = 10^ 
cm“^ for the sun) but much smaller than the photo¬ 
sphere. Therefore the thickness of the emitting region 
should be not higher than the corona of the companion. 
On the other hand the radiation length of the relativistic 
electrons is quite long: 716Af {Z{Z + 1) ln(278/ \/\Z)) ^ 
4 X 10^^ cm for the coronal gas of np ^ 10^cm“^. There¬ 
fore these electrons can radiate synchrotron X-rays even 
in the immediate area of the photosphere. 

If the hard X-rays originated in the shocked pulsar 
wind on the surface of the companion star as is discussed, 
the composition of the pulsar wind can be constrained 
based on the downstream flow velocity. Gonsidering 
the Rankin-Hugoniot relation for magnetized electron- 
positron plasma for small a, we can find the downstream 
flow velocity as a function of the magnetization parame- 
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ter (j (iKennel fc CQrQiiitilll984j) . 


where U 2 = ^ 2 ^ 2 - Finally we find the downstream flow 
velocity, (32 = ^(1 + 9<j)/(7 — 9cr) in units of the light 
speed. Note that the equation predicts the lower limit of 
the shocked pulsar wind, ^^2 = l/\/7 ^ 0.38 at cr = 0. On 
the other hand, the flow velocity of P 2 = 0.43 ±0.15 from 
the X-ray light curve indicates that the pulsar wind is al¬ 
ready in the particle dominant state, i.e., a ^ 0.028, at 
the surface of the companion, only 1.1 x 10^^ cm apart 
from the pulsar. Note that the post-shock flow veloc¬ 
ity has an inverse-correlation with the inclination angle 
that affects the resultant magnetization parameter. If 
we adopt an error range of inclination angle, 50° ^ 60°, 
the magnetization parameter can range from 0.008 (for 
i = 60°) to 0.1 (for i = 50°). The upper limit of cr < 0.1 
is still much smaller than that calculate d for the similar 
system of PSR J1023±0038 reported bv iBogdanov et all 
mnh , in which a thick post-shock emitting region was 
assumed. Past theoretical studies claimed that the pulsar 
wind is highly magnetized at the beginning just around 
the light cylinder. However, the observed cr is far smaller 
than 1, and t his discrepancy has been called the cr para¬ 
dox. Recently I Aharonian et al reported evidence 

of particle acceleration just around the pulsar based on 
gamma-ray pulse analysis of the Crab pulsar. The above 
argument of a in this work is consistent with their result 
and may be one of the cr parameters measured at the 
nearest distance from a pulsar. Furthermore this may 
imply that the old MSPs and the young radio pulsars 
have the same particle acceleration mechanism. 

One of the intriguing features of this object is the 
double-peak structure in the 2—4 keV X-ray light curve. 
Since the dip is observed only in the softer band (Fig¬ 
ure ED, this feature can be interpreted as absorption of 
soft X-rays. If we accept that the non-thermal X-rays 
originate in the shocked pulsar wind on the surface of 
the companion as discussed above, the absorber must be 
lying around the pulsar, because the dip appears at the 
optical maximum phase (= the superior conjunction). 
In addition, the orbital inclination angle i ^ 55° requires 
that the absorber must be extended far from the orbital 
plane (^ 8 x 10^^ cm) . 

As iRnv et aP (|2Q14l ) claimed, mass accretion on the 
pulsar seems to continue intermittently, however the ar¬ 
gument of the propeller effect rules out a powerful ac¬ 
cretion. This is consistent with the lack of a bright disc 
component in X-ray spectra. In addition to the pro¬ 
peller effects, the measured companion radius is about 
1/2 of the Roche lobe, therefore the stellar material can¬ 
not overcome the Lagrange point LI to accrete on the 
neutron star. To feed mass to the pulsar, explosive ejec¬ 
tion faster than ^ 1000 km s“^ may be required for the 
mass to escape from the gravity potential. 

In addition, as inferred from the discussion on the heat¬ 
ing efficiency, the pulsar wind may be concentrated on 
the orbital plane, which would prevent the accretion flow. 
If the stellar material that is blown away by the equa¬ 
torial pulsar wind is drifting around the pulsar’s polar 
regions, and it may accrete from the poloidal direction. 
On the other hand, the running periodogram shown in 


Figure [TOl indicates that the dip structure related to the 
second harmonic component is unstable, and appears to 
have disappeared ^ 120 ks from the start time of the ob¬ 
servation. To confirm the scenario in which the stellar gas 
is drifting around the pulsar and accretes intermittently, 
simultaneous X-ray and radio observations are required. 

5. CONCLUSION 

We presented multi-wavelength observations of a newly 
found black-widow binary system 2FGL J2339.6-0532 
covering near-infrared to X-ray regimes. Thanks to a 
wide wave waveband and long coverage, we successfully 
obtained a phase-resolved SED that enabled us to con¬ 
strain the orbital parameters more precisely. The ob¬ 
tained SED seems consistent with past studies, and we 
calculated an inclination angle of i ^ 55°, while tak¬ 
ing into account the results of recent radio observations. 
Based on the argument of the propeller-effect argument, 
the lower limit of the acrretion rate was estimated to 
be m > 1.3 X 10^^ g which is about five orders of 
magnitude larger than that expected from the observed 
X-ray luminosity. In addition, the estimated orbital pa¬ 
rameters imply that the companion’s radius is only 1/2 
of the Roche lobe, making it difficult to feed the pul¬ 
sar continuously. We also obtained an Ra image of the 
vicinity of the target and could not detect any diffuse 
structure with a 3cr detection limit of < 8.7 X IQ-i^ erg 
s“^ cm“^ arcsec“^. We therefore conclude that the tar¬ 
get does not have a bow-shock nebula brighter than the 
Ra nebula around PSR B1957±20. 

In the X-ray regime, we discovered a steady, soft X- 
ray component below 1 keV, which seemed to originate 
from the neutron-star surface, and which showed no ev¬ 
idence of an accretion disc. On the other hand, the 
hard component above 2 keV showed periodic modula¬ 
tion synchronized with the orbital motion, implying that 
the hard X-rays originate in the shocked pulsar wind 
near the companion surface. The observed X-ray lumi¬ 
nosity is comparable to the bolometric luminosity of the 
heated hemisphere of the companion. This means that 
the heating efficiency should be smaller than 0.5 and the 
pulsar-wind distribution should be anisotropic. Adopt¬ 
ing the above scenario we estimated the magnetization 
parameter of the pulsar wind a ^ 0.03, based on the 
Rankin—Hugoniot relation. This implies that the pulsar 
wind is already in the particle-dominant state at a dis¬ 
tance of 1.1 xlO^^ cm from the pulsar. Moreover, we 
also investigated the time variability of the modulation 
pattern by using the running periodogram of X-ray light 
curve, wherein we detected a weakening of the modu¬ 
lation pattern. This may be related to an intermittent 
weak accretion or an unstable pulsar-wind activity, both 
of which can cause temporal extinction of the radio-pulse 
emission. 
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Table 4 

Summary of optical observations with OISTER 


Telescope (Instrument) 

Diameter 



September 2011 





October 2011 

22 

23 

24 

27 

28 

29 

30 

1 

4 6 7 

Pirka telescope (MSI)®' 

160cm 


BVRit 

BVRit 

V 

V 





Gunma Astronomical Observatory (GIRCS)^ 

150cm 





J 





MITSuME-Akeno telescope (Tricolor Camera)*^ 

50cm 

gRI 

gRI 

gRI 

gRI 

gRI 

gRI 

gRI 



Kiso Schmidt (2kCCD)‘^ 

105cm 






RHo; 




Kyoto Sangyo University 1.3m telescope (CCD)® 

130cm 





gz 





Nayuta telescope (NIC)^ 

200cm 

(JHKg) 









MITSuME-OAO telescope(Tricolor Camera)® 

50cm 


gRI 


gRI 

gRI 



gRI 

gRI gRI 

OAO 188cm telescope (ISLE)® 

188cm 




J 






Bisei Spaceguard Center Im telescope (Optical CCD)^ 

100cm 




(r) 

(r) 





Kanata telescope (HOWPol)j 

150cm 




R 

R 




V B 

Kagoshima University Im telescope (IR CCD)J 

100cm 




(J) 

(J) 



(J) 


Murikabushi telescope (Tricolor Camera)® 

105cm 




gRI 


gRI 

gRI 



Kottamia Astronomical Observatory (Optical CCD)^ 

188cm 




R 

R 

R 

B 



IRSF 1.4m telescope (SIRIUS)^ 

140cm 




JHKg 

JHKg 

JHKg 



miniTAO (ANIR)^ 

104cm 






RJ 

BJ 




Note. — The letters on the table describe the observed energy bands. Lower case and upper case characters correspond to the SDSS system and the 
Johnson —Cousins system, respectively. The observations shown in parentheses did not produce usable data due to the weather condition. 


t 

a 

b 

c 

d 

e 

f 

g 

h 


Calibration observation for the field photometry. 
Nayoro, Hokkaido Pref., Japan IWatanabe et al.ll201~^ 


Takayama, Gunma Pref., Japan 
Hokuto, Yamanashi Pref., Japan 


ITakahashi et al 

iKotani et ^ 




Kiso, Nagano Pref., Japan mK et al DUMH) 

Kyoyo, Kyoto Pref., Japan I http://www.kyoto-su.ac. jp/kao/I 
Sayo, Hyogo Pref., Japan I http://www.nhao. jp/en/I 
Asaguchi, Okayama Pref., Japan lYanagisawa et al.ir2UlUI ') 


Asaguchi, Okayama Pref., Japan i 


^ Bisei, Okayama Pref., Japan I http://www. spaceguard.or. jp/bsgc_jsf/pamphlet/index.htm I 


J Higashi-hiroshima, Hiroshima Pref., Japan JKawabata^^t_aJj[200^ 

Satsuma-sendai, Kagoshima Pref., Japan |http://milkyway.sci.kagoshima-u.ac.jp/lm/sys/index.html | 
^ Ishigaki, Okinawa Pref., Japan iKotani et an i200^') 

^ Kottamia, Egypt lAzzam et al.||2014|) 

™ Sutherland, South Africa INagavama et al.|[^U3| ') 

Atacama, Chile onis Motohara et aT1l2Ul(l') 



























